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ABSTRACT 

We present a detailed investigation into which properties of cold dark matter halos make them 
effective strong gravitational lenses. The cross sections for giant arc formation of 878 clusters from 
a high-resolution N-body simulation of the ACDM cosmology are measured by ray tracing through 
13,594 unique projections. We measure concentrations, axis ratios, orientations, and the amount of 
substructure of each cluster, and compare the lensing weighted distribution of each quantity to that 
of the cluster population as a whole. We find that NFW profiles provide just as good a fit to lensing 
clusters as they do to the total cluster population; however, the concentrations of lensing clusters are 
on average 34% larger than the typical cluster in the Universe. Despite this bias, the anomalously 
high concentrations of lensing clus ters (c > 14) recently measured by several groups from combine d 
strong and weak lensing analyses ijKneib et al J 120031 iGavazzi et al.ll2003t iBroadhurst et al.ll2005bj) . 
appear to be inconsistent with the concentration distribution in our simulations, which predict < 2% 
of lensing clusters should have concentrations this high. No correlation is found between strong lensing 
cross section and the amount of substructure, indicating that the population of cluster lenses is no 
more relaxed or disturbed than typical clusters in the Universe. Lensing clusters tend to have their 
principal axis aligned with the the line of sight: the median angle is |cos6*| = 0.67. We introduce 
several different types of simplified dark matter halos, and use them to isolate which properties of 
CDM clusters make them effective gravitational lenses. Projections of halo substructure onto small 
radii and the large scale mass distribution of clusters do not significantly influence strong lensing 
cross sections. The abundance of giant arcs is primarily determined by the mass distribution within 
an average overdensity of ~ 10, 000. The clumpy cores of dark matter halos result in ~ 25 — 60% more 
giant arcs than smooth ellipsoids of the same total mass. A multiple lens plane ray tracing algorithm 
is used to show that projections of large scale structure increase strong lensing cross sections by a 
modest amount < 7%. We revisit the question of whether t here is an excess of g iant arcs detected 
for high redshift clusters in the Red-Sequence Cluster Survey ijGladders et al.ll2003|) and find that the 
number of high redshift (z > 0.6) lensing clusters is in good agreement with ACDM, although our 
simulations predict more low redshift (z < 0.6) lensing clusters than were observed. 
Subject headings: dark matter - galaxies: clusters: cosmology: theory - methods: numerical - clusters: 
general - large scale structure of the universe - gravitational lensing 



1. INTRODUCTION 

Strong gravitational lensing by clusters of galax- 
ies provides a unique laboratory for studying the 
small scale dark matter distribution of the largest 
collapsed structures in the Universe. Because strong 
lensing directly probes the gravitational potential, 
it is free from the assumptions which plague other 
techniques which probe the mass distribution on 
compara ble scales. For examp l e, dynamica l mea- 
sures l|Nataraian fc Kneibl Il996t iKelson et, all 120021) 
assume that dark halos are virialized systems and 
models of the X-ray te mperature profile of the intr- 
acluster medium (e.g. lArabadiis. Bautz. fc Garmirel 
20021 lEttori. Fabian Allen, fc Johnstone! I2002t 
Lewis. Buote. fc Stockd l2003j) typically assume hy- 
drostatic equilibrium. 

Up coming X-ray llRomer et all l200"ll lEbeling et al. 
1 20011) Sunyaev-Z eldovi ch iCarlstrom Holder, fc Reese 
120021 iKosowskvl F2003t iSchwan et all 12003ft . and opti- 
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cal (jGladders fc Yeel l2004t) cluster searches will dra- 
matically increase the number of clusters of galax- 
ies known. Deep follow up imaging of these clus- 
ters will discover hundreds of new giant arcs. Such 
statistical samples of arcs will allow measurement of 
the giant arc abundance as a function of cluster red- 
shift, the radial distributions of arcs from the clus- 
ter center, and the distribution of relative angles be- 
tween arcs in multi-arc systems. Comparisons of these 
quantities to expectations from ray tracing simulations 
of clusters i n a ACDM cosmology have already been 
carried out llBar telmann et al.lll 998t iWambsganss et al. 
| 2004a| iDalal Holder, fc Hennawil 12004: Ho fc White 
l2004HLi e t^jdj[2005'l using small samples of ~ 20 known 
arcs fsee lSand et al.ll20051 for a recent compilation) . 

Armed with future statistical samples, we will be able 
to measure the distributions of cluster lens properties, 
in addition to the simple 'one point' statistics of the 
arcs. Detailed modeling of image positions in clus- 
ter lenses can determine a set of best fit parameters 
which describe the distribution of dark matter in each 
cluster lens (e.g.lTvson. Kochans ki. fc dell'Antonioll998t 
Smit h et all 120011 ISand et alJ 120041 IBroadhurst et alJ 
l2005ajr ~and even stronger constraints can be obtained 
when strong lensing information is combined with larger 
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scale weak lensing measurements (Kneib et al.l l2003t 
lOavazzi et al.ll2003t iBroadhurst et al.ll2005blL These pa- 
rameters might include the slope of the mass profile, its 
concentration, or the projected ellipticity of the cluster. 
A comparison of the observed distribution of these pa- 
rameters to the expected structures from cosmological 
simulations of a ACDM Universe will provide a strict 
test of our current theory of dark matter and structure 
formation. 

In this context, the questions naturally arise: Do 
strong lensing selected clusters constitute a fair sample of 
clusters in the Universe? Do biases exist with respect to 
halo concentration, triaxiality, preferential alignment, or 
the amount mass in substructure? Clearly, these biases 
must be taken into account before a comparison can be 
made between the observed distribution of cluster lens 
properties and the properties of numerically simulated 
clusters. Additionally, by comparing the distribution of 
lensing cluster properties to the cluster population as a 
whole, we can hope to gain insights and intuition about 
what makes a CDM cluster an effective strong lens. 

In this paper we statistically characterize the proper- 
ties of the cluster lens population. Using the ray trac- 
ing technique discussed lDalal. Holder, fc HennawJ 1)2004 
henceforth DHH), we compute the lensing cross sections 
of a large sample of clusters from a cosmological N-body 
simulation. We measure the properties of each cluster 
and compare the lensing weighted distribution of clus- 
ter properties to the cluster population as a whole. We 
also attempt to isolate which properties of CDM clusters 
make them effective gravitational lenses. In this context 
we introduce 'Analog Halos,' which are simplified dark 
matter halos which retain one or more of the proper- 
ties of the original simulated clusters. As an aside, we 
also quantify the effect of projections of large scale struc- 
ture between the source plane and the observer on the 
abundance of giant arcs. Armed with lensing cross sec- 
tions for clusters in a large simulation volume, we revisit 
the question of whether there is an excess of giant arcs 
detected for high redsh ift clusters in the Re d-Sequence 
Cluster Survey (RCS) ^ladders et al]l200lh . We were 
unable to answer this question definitively in DHH be- 
cause the cosmological volume simulated in that study 
was too small (see also FWambsg anss et al.ll20"04aD . 

The cosmological N-body simulation and our method- 
ology for measuring the cluster parameters and quanti- 
fying substructure are described in § We review the 
ray tracing algorithm in § and investigate the effect of 
additional matter in the light cone in §0J In §0 we intro- 
duce the 'Analog Halos' and examine which properties of 
CDM halos make them effective gravitational lenses. We 
compare the statistical properties of cluster lenses to the 
total cluster population in § El • The abundance of giant 
arcs in the RCS is revisited in § Q and we summarize 
and conclude in § IH1 

2. SIMULATED CLUSTERS 

In this section, we describe the simulated clusters 
which are the input to our ray tracing algorithm. For the 
analysis in §[3]and §|S]we will require various properties of 
clusters, such as triaxiality, the amount of substructure, 
and cluster concentration both in three dimensions and 
for two dimensional projections. First we describe the 
N-body simulations and how the clusters were identified, 



then we describe how cluster properties were measured. 

2.1. N-body simulations 

In order to simulate the strong lensing effects of a Uni- 
verse filled with dark matter, we used clusters drawn 
from a cosmological N-body simulation. The simulation 
w as performed with the T ree-Particle-Mesh (TPM) code 
of lBode fc Ostrikerl <|2003j) . TPM uses the Particle-Mesh 
method for long-range forces and a tree code for sub- 
grid resolution; individual isolated, overdense regions 
are each treated as a separate tree, thus ensuring effi- 
cient use of parallel computers. The following cosmo- 
logical parameters were used: matter content S1m 
0.3. cosmological constant = 0.7, Hubble constant 
Hq = 70 km s _1 Mpc -1 , linear amplitude of mass fluc- 
tuations (T8 = 0-95, and primordial power spectral index 
n s =l. These parameters are consistent (wi thin la) of the 
WMA P derived cosmological parameters l)Spergel et alJ 
2003). The simulation volume is a periodic cube with 
comoving side length of L = 320 hr 1 Mpc containing 
N = 1024 3 particles, so the particle mass is m p = 2.54 x 
10 9 /i _1 Mq. The cubic spline softening length was set to 
e = 3.2/i _1 kpc. The small softening length, small parti- 
cle mass, and large volume make this simulation ideal for 
studying strong lensing by clusters; this simulation was 
previously analyzed bv lWambsganss et alJ l)2004albfl . We 
used outputs at seven different 'snapshots', spaced a co- 
moving length L apart and covering the range of redshifts 
over which the critical density is low enough to produce 
an appreciable amount of strong lensing. These red- 
shifts are z = 0.17,0.29,0.41,0.55,0.70,0.87,1.05,1.26, 
and 1.49. 

A "friends-of- friends" (FOF) group finder 
l)Davis et, alJ E?.85) with the canonical linking length 
of 6 = 0.2 was applied to each particle distribution to 
identify cluster size dark matter halos. For each cluster 
with a FOF group mass above Mfof > 10 14 /i _1 M , 
all the particles within a 5 hr 1 Mpc sphere about the 
center of mass were dumped to separate files and used 
as inputs to our ray tracing code. 

2.2. NFW Profile Parameters 

For each cluster identified by FOF, the gravitational 
potential was computed and the center of the cluster was 
defined to be the most bound particle in the 5 h^ 1 Mpc 
particle dump. The radial mass profile of the cluster 
about this center was calculated, and the virial radius 
was defined as the radius of the innermost particle at 
which the average density interior to it was greater than 
or equal to A v j r (z) times the cosmic mean matter den- 
sity where the virial overdensity is A v ; r w 187r + ®+ x x ~ 39x 
l)Brvan fc Normanlll998j) . with x = Q m (z) - 1. 

We computed the number of particles in 50 logarithmi- 
cally spaced radial bins between the 100th particle from 
the cluster center and the virial radius. The smallest 
radial bin was placed at the position of the 100th par- 
ticle because this roughly defines the scale at which the 
mass d istribution is unaff ected by two body relaxation 
effects ijPower et al.ll2003l . This scale exceeds the soft- 
ening length for all of the clusters we considered and 
provides a stricter limit on the scales we resolved. 

We then fit the number of particles in each bin d/Vj to 
the number profile corresponding to the universal density 
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Fig. 1. — Example of a cluster with a large amount of substructure inside the virial radius. Left: Three dimensional density distribution 
of the cluster. Surface density projections on each coordinate plane are also shown to aid visualization. The spatial axes are in units of 
comoving h _1 Mpc and both density (2-d and 3-d) color maps are logarithmically stretched. Right: Radial density profile (points) and best 
fit NFW profile fit (solid green line), which has a reduced \ 2 = f.28. The y-axis has units of h _1 Mq/pc 3 . The inner and outer dashed 
vertical lines indicate the locations of the scale radius and virial radius for the NFW fit, respectively. The lower panel plots the ratio of 
Ap, which is the difference between the profile and the fit, to a p , the bootstrapped density variance, for each radial bin. The bootstrapped 
variances properly down weight the bump in the outer radial bins caused by the massive substructures. Had we used the Poisson errors 
alone the radial bump would have had a catastrophic effect on the fit because 1/y/N is very small for these outer bins. 



profile of iNavarro. Frenk. fc White! l)1997|) (henceforth 
NFW). The number of particles in each bin dNi ~ r 2 pdr 
is much less noisy than the density and furthermore, it 
is dNi which should obey Poisson statistics if counting 
errors are the only source of noise in each measurement, 
allowing us to assign sensible errors (see discussion be- 
low). 

Thus we fit the quantity 

dNi(r s , M c ) = M c [A( Xi ) ~ A{x^ x )\ (1) 
for the two parameters Mc and r s , where 

A(x)=Hl + x)--^-. (2) 
1 + x 

and x = r/r s . In terms of the familiar NFW parameters 
Mc = 47rp s rf/ra p , but we opt to fit for Mc because the 
profile is linear in this quantity. 

The parameters are determined by minimizing the re- 
duced x 2 penalty function 

2 1 ^ [dN-dN^Mc f 

z— 1 ' 

Our computation of the variance for each logarithmic 
bin, of, warrants further discussion. The deviations be- 
tween our simulated clusters and an NFW profile can 
be caused by both counting errors, described by Poisson 
statistics, and massive substructures or asymmetries in 
the cluster. Without taking the latter source of error 
into account, a bump in the radial profile of a halo will 
have a catastrophic effect on the NFW fit if it is in the 
outer region of the halo, because the formal Poisson er- 
rors for these radial bins will be very small. Previous 



workers have dealt with this prob lem by removing th ese 
substructures from the halo (e.g. Uing Sz Sutdl2002|) or 
truncating the profile and the fit wherever these bumps 
occur (e.g. IBullock et al1l200H) . However, in the context 
of strong lensing, neither of these approaches is appro- 
priate. We wish to measure a three dimensional profile 
which is as closely related as possible to the two dimen- 
sional surface density profile probed by strong lensing 
observations. Strong lensing is sensitive to the total pro- 
jected mass, and since it is not possible to remove or 
truncate substructures observationally we opt not to do 
so in the simulations either. 

Instead, we define a modified variance which incorpo- 
rates halo asymmetry but which reduces to the Poisson 
variance for a spherically symmetric halo. Each radial 
bin is divided into 12 equal area sectors using Healpix 5 , 
and we estimate the variance via a 'bootstrap' 

^^-"fk-^y, ( 4) 

Tlpix j—i V n pixJ 

where n p ; x = 12 is the number of Healpix sectors and 
driij is the number of particles in the ith radial bin which 
land in angular sector j. 

The right panel of Figure ^ compares the measured 
density profile to the best fit NFW profile 6 , for a cluster 
from our simulation which is depicted in the left panel. 
The massive substructure interior to the virial radius 

5 http:/ /www. eso.org/science/healpix 

e Although we fit the number profiles using dNi and <7j, we plot 
density profiles p(r) and corresponding propagated errors a p for 
the sake of visualization. 
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causes a bump in the radial density profile, however, the 
bootstrapped errors correctly down weight these points 
and minimize their influence on the fit. 

The distribution of the reduced x 2 for the 878 clusters 
in our fiducial snapshot is displayed in Figure The 
fact that median value of this distribution is near unity 
gives us confidence that the NFW profile is providing a 
good fit to our simulated clusters and our bootstrapped 
(7j are accurately describing the Poisson errors as well as 
the deviations from symmetry. 

Once Mq and r s are determined, all the NFW param- 



eters r s , p s 



and M v ir can be determined from 



simple algebraic relations. Although the NFW fit deter- 
mines the virial mass of the cluster, we instead adopt the 
non-parametric virial mass given by the mass interior to 
the virial radius defined by the overdensity criterion A v i r . 

We also perform fits of the surface mass density to 
a two dimensional projected NFW profile, for various 
projections through the mass distribution of each clus- 
ter. These two dimensional fits are of interest because 
the projected mass distribution is the quantity acces- 
sible to observations of strong lensing. We repeat the 
fitting procedure described above for surface mass den- 
sity by computing the azimuthally averaged radial pro- 
file in 50 logarithmically spaced bins between 0.02 and 
2 h^ 1 Mpc. The same bootstrapping procedure (eqn. 
is used to compute the errors, where we use 20 az- 
imuthal bins to estimate the error for each radial bin. 
We minimize the x 2 of a fit to an analytical expres- 
sion for the projected N FW profile l(Bartelmannl Il99fj 
iWright k Brainerdl2000|) . Note that this formula for the 
projected NFW profile assumes the profile extends to in- 
finity, which results in a slight bias of 6% (see § I6.3JI in 
the concentrations measured. For the sake of simplicity, 
we use this analytical form rather than fitting to the ex- 
act numerical expression for a projected truncated NFW. 
This small bias does not have a significant effect on our 
conclusions. 

2.3. Substructure 

Although complicated algorithms exist to identify sub- 
struc t ures in N-body simulati o ns Csee e.g.lKravtsov et al. 
20031 iDe Lucia et all 12004) iWeller. Ostriker. k Bode 



2004; IGill. Knebe. fc CxibsoJ^OOl " for our purposes we 
are only interested in identifying substructures massive 
enough to influence the strong lensing cross sections of 
massive clusters, for which even a very simple procedure 
will suffice. In addition, there is no need to use veloc- 
ity information to unbind particles as is done in other 
algorithms since the energy of the matter distribution is 
irrelevant to its lensing effect. To this end, we identify 
substructures by rerunning the FOF algorithm on the 
particles in each cluster dump, but with a smaller link- 
ing linking length of b = 0.05. This causes a cluster halo 
to fragment into one or more dense subhalos. The most 
massive subhalo will be centered on the potential mini- 
mum or center of the cluster, and we designate this as 
the cluster core. We define as a substructure any addi- 
tional subhalos (besides the most massive) identified by 
the FOF algorithm with a mass M > 10 12 h' 1 M©. This 
is safely larger than the smallest subhalo resolved by our 
particle mass ~ 100m p = 2.5 x 10 11 h^ 1 M . 

A linking length b = 0.05 groups particles within an 
isodensity contour of ~ 4, 000 times the background den- 



sity l|Lacev k Cold Il994|) . For an NFW cluster, the 
smaller linking length truncates the cluster at a radius 
~ r io,ooOi where the average overdensity is ~ 10,000 
times the mean density. For, M v - U ~ 5 x 10 14 h~ l M , 
Cyir = 6, and z = 0.4, rio,ooo ~ 300 h" 1 kpc (comoving), 
which projects to an angular separation of ~ 1'. 

We define three different statistics which quantify the 
amount of substructure in a cluster. Two of the three 
can be written as 
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N sub 



^norm n 



(5) 



where Mi refers to the mass of the i subhalo, M noTm 
is a normalizing mass, and the sum extends over all sub- 
structures enclosed within the virial radius of the halo. 
Note that the sum excludes subhalo i = 1, which as noted 
above corresponds to the cluster core. 

We explore two different normalization masses 
M norm = M C ore and M n0 rm = Afvir, corresponding to the 



two statistics M su b- 



and M, 



sub— vir j 



respectively. The 



statistic M su b-corc is normalized by the mass of the dense 
core which is the mass relevant to lensing. This could be 
biased low for lensing clusters since they might have more 
massive cores which is why we also consider M su b- V ir- 
Our third substructure statistic is M su b-2 = M^jM^, 
where M2 is the second most massive sub-halo. This 
statistic provides a measure of the binarity of the clus- 
ter, and is less sensitive to lower mass halos which may 
dominate the total mass in substructure because they are 
more abundant. 

2.4. Triaxiality 

Previous work has shown that for CDM halos, 
the degree of triaxiality increase s (i.e., the axis ra- 
tios decrease) toward the center (IBa rn^sfcl^fstatdiiai 



T^llWarren et al.ll992H.Iing k Sutoll2002HSchulz eTaf 



2005). For our purposes, we arc interested in the triaxi- 
ality of the mass within the dense core which is primarily 
responsible for the strong lensing. As mentioned in the 
previous section, we identify the most massive subhalo 
selected by the FOF algorithm using a smaller linking 
length with the cluster core. We compute the principal 
axes and axis ratios by diagonalizing the normalized mo- 
ment of inertia tensor for the particles which are linked 
into this core subhalo 



Q 



1 



E 



(G) 



where N COTC is the number of particles grouped into the 



core subhalo. If 



> b 2 > 



are the principal com- 



ponents of Q, we compute the axis ratios 172 = b/a and 
q 3 = c/a for each cluster. The eigenvectors of Q give 
the orientation of the principal axes V a , Vb, V c . Note 
that by restricting attention only to particles which lie 
in the FOF group, we are computing the axis ratios of 
the particles within an isodensity contour, similar to the 
approach of Ding k Sutol ((2002). Other workers have 
computed the axis ratios of particles within an ellipsoid 
ijBarnes k EfstathioullT987t IWarren et al.lll992|) defined 
by an iterative procedu re, but these iterations are not 
guara nteed to converge l)Jing k Sutcll2002l ISchulz et alJ 
12005(1 which is why opt for the method used here. 
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3. STRONG LENSING CALCULATION 

We compute the strong lensing properties of our simu- 
lated clusters using ray-tracing; cross sections and arc 
statistics are determined via Monte Carlo integration 
over the background source plane. We refer the reader 
to DHH for details of the ray tracing and the Monte 
Carlo calculation. The ray trace and Monte Carlo are 
repeated for each projection through the cluster and for 
each source plane considered. We ray trace five source 
planes at z s — 1.0, 1.5, 2.0, 3.0, and 4.0. The num- 
ber counts of background galaxies are collapsed into 
five source redshift bins centered on these source planes 
given by the intervals [0.75,1.25], [1.25,1.75], [1.75,2.5], 
[2.5,3.5], and [3.5,6.0], respectively. Because the critical 
density for strong lensing is a slowly varying function of 
source redshift, this binning should not introduce signif- 
icant errors in our calculation. Given the cross sections 
for each source plane and the number counts of back- 
ground galaxies in each bin, the number of giant arcs 
behind each cluster orientation can be computed as a 
function of li miting magn i tude o r surface brightness. 

DHH and IHo fc W hite ( 2004) found that the lensing 
cross section for a single cluster varies enormously as a 
function of orientation (see Figure 3 of DHH). A fac- 
tor of ~ 20 spread in lensing cross section was not un- 
common and the histogram of cross sections for differ- 
ent orientations was quite skewed. To compute reliable 
mean cross sections for a cluster, we must thus aver- 
age over many projections to appropriately sample the 
distribution of cross sections. This orientation averag- 
ing is most important for the rarest, most massive clus- 
ters in the simulation volume that dominate the total 
lensing optical depth. We computed mean lensing cross 
sections by averaging over 125 orientations for all clus- 
ters with Mfof > 10 147 hr 1 Mo, 31 orientations for 
clusters in the mass interval 10 143 h^ 1 M < Mfof < 
lO 14 7 Yr 1 Mq, and 3 orientations for all clusters in the 
range 10 14 hr 1 M Q < M FO f < 10 14 3 hr 1 M . 

Given the cross sections for each cluster orientation, 
the total number of arcs caused by the clusters in the 
snapshot at redshift Zj 



\ ^ 



L 3 



<Tijk(> 0). 



(7) 



Here Oijfc(> 0) is the orientation averaged cross section 
for forming giant arcs with angular separation larger than 
9 for the ith cluster in snapshot Zj with source plane z k \ 
f2 is the solid angle of the survey under consideration, 
L 3 is simulation volume, and n k is the number density 
of background galaxies in the bin about source plane z k . 
The quantity Vj is the effective volume represented by 
the snapshot 



dV 
dfldz 



dz, 



(8) 



where is the cosmological volume element. Hence, 
we take the clusters in the snapshot at redshift Zj to be 
representative of the lensing rate in the volume of the 
light cone over the redshift interval [ Zj+ ^ J ~ 1 , Zj+ * J+1 ], 
Finally, in what follows we will often use the effective 



cross section of a cluster, which we define as 

Efe n fcCT fe (> 6) 



o- e ff(> <?) 



Efc n k 



(9) 



We defer a discussion of the counts of background 
galaxies in each source redshift bin, n k , to § 17.11 

3.1. Adding Brightest Cluster Galaxies 

Several groups have considered the degree to which 
cluste r galaxies alter the cross s ection for producing giant 
arcs. IMeneghetti et all l|2000fl and iFlores et all ( 2000) 
studied the influence of cluster galaxies (besides the 
BCG), and found them to be generally unimportant. The 
effects of the cen tral brightest cluster galaxy (BCG) were 
conside red by IMeneghe tti. Bar telmann. fc Mosc ardini 
( 2003b), DHH, and IHo fc White! l)2004|) . These stud- 
ies found that BCGs increase the lensing cross section 
for arc separations small enough that the mass enclosed 
within the arc radii ha s a significant baryq nic component. 
Specifically, DHH and lHo fc Whitel pool found that for 
arc radii > 10", the difference in the number of arcs for 
pure dark matter clusters versus clusters including BCGs 
was not dramatic. 

In order to account for the effects of BCGs on the 
giant arc cross sections, we artificially add baryons to 
the centers of each cluster by 'painting' BCGs onto the 
dark matter surface density. Similar to the procedure in 
DHH, we model the mass profile with singular isothermal 
spheres (SIS), which accurately represent strong lens - 
ing by elliptical galaxies (see e.g. iKochanek et afl feOQO'). 
There it was found that the BCG enhancement of the 
cross section depended sensitively on the concentration 
of the central galaxy, parameterized by the velocity dis- 
persion of the SIS; whereas varying the total mass of the 
central galaxy (i.e. the extent of the SIS) had a negligible 
effect. We thus choose to keep the mass of the central 
galaxy to be a fixed fraction Mbaryon = 0.003 A/fof of 
the mass of the dark halo. 

In the absence of a theory of BCG formation, we use a 
simple prescription for assigning velocity dispersions to 
the B CGs, scaling it with th e mass of the dark matter 
halo. lEdge fc Stewartl l)1991fl found a strong correlation 
between the optical luminosity of BCGs and the X-ray 
temperature of a large sample of galaxy clusters, with 

scaling £bcg °c T^ 10 . Because BC Gs lie on the fun- 
dame ntal plane of elli ptical galaxies llOegerle fc Hoessel 
Il991|) . they obey the iFaber fc Jackson! l)1976j) relation 
L oc it 4 . Combining these scaling relations with the 
mass-temperature relation for X-ray clusters, Tx 
M 2 / 3 , gives a oc M 2 / 15 . Normalizing this relationship 
to the Coma cluster gives 



300 km s" 



Mi 



FOF 



lO 15 /^ 1 M 
velocity 



2/15 



(10) 



where we took the velocity dispersion of 
the BCG in Coma to b e a = 323 km s _1 
ijFisher. Illingworth. fc Franxl |^995i) and converted 
the X-ray temperatu re of Coma, Tx = 8.0 keV 
(|Edge fc Stewartlll991|) . to a FOF gr oup mass (b = 0.2) 
usin g the prescripti on described in iWhitd l)2001|) (see 
also iHu fc Kravtsovll2lj03l) . 
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Fig. 2. — Cumulative distribution of giant arc separations for 
both single plane (dashed red) and multi plane (solid black) ray 
tracing simulations. Our full ray-tracing simulations show that line 
of sight projections have a modest effect, < 7%, on the statistics 
of giant arcs with separations 9 > 10". 



4. QUANTIFYING THE EFFECT OF MATTER IN THE 
LIGHT CONE 

A question of recent interest has been the importance 
of projections of large- scale structure in th e line of sight 
towards strong lenses. IWambsg anss et all l|2004b|) have 
suggested that projections of large scale structure signifi- 
cantly (~ 20—30%) increase strong lensing cross sections, 
using heuristic arguments. In this section, we quantify 
the effect of projections on giant arc cross sections using 
full ray tracing simulations. 

We generalized our ray tracing code to compute lens- 
ing by matter at different reds hifts using the multi- 
ple lens plane algorithm (see e.g. iSchneider et al.l fl992'). 
This allows us to include the strong lensing effect of 
the large scale structure in the light cone between the 
source plane and the obs e rver. This code was used in 
iDalal. Hennawi. &: Bodel (|2005) to study the effect of 
projections of large scale structure on strong lensing cos- 
mography, and a full description can be found there. 
Here, briefly summarize the essential features of the code. 

We sliced through our 320 h^ 1 Mpc simulation cube to 
produce lens planes spaced every 160/i _1 Mpc with an- 
gular size roughly 20', so that the physical side length 
of the projected planes telescopes with increasing dis- 
tance from the observer at z = 0. A total of 243 pairs 
of planes were produced for each of 19 redshifts between 
z = — 6.37. The cosmic mean density was subtracted 
from each plane, and the lightcone was tiled by randomly 
selecting a pair of planes for each redshift, resulting in 
a total of 38 lens planes per realization. This procedure 
and the simulation outputs a re identical to those used by 
IWambsganss et all l)2004alblh 

Ideally, we would compute arc cross sections simply by 
generating many such realizations and ray tracing over 
the entire 20' field. However, because the computation 
of reliable mean cross sections requires averaging over 
a large number of different cluster orientations, this ap- 
proach would not be computationally feasible: we would 
spend the vast majority of our time ray tracing blank 



regions of sky far from any clusters which would contain 
no giant arcs. To circumvent this practical constraint, we 
simply insert an extra lens plane containing the massive 
clusters for the lens redshift or snapshot that we consider. 
We ray trace only a 4 /i -1 Mpc FOV around each cluster 
and project through the same number of different orien- 
tations as described in the previous section. However, 
each orientation now also includes a different random re- 
alization of the 38 lens planes to account for the large 
scale structure along the line of sight. This allows us 
to compute the total orientation-averaged strong lensing 
cross section for the clusters in this snapshot, but with 
the additional effect of large scale structure in the light 
cone. A direct comparison can then be made to the ray 
traces which included the clusters alone, as the orienta- 
tion averaging is identical for both cases. We will refer to 
the ray traces which include large scale structure along 
the line of sight as the 'Multi-Plane' simulation and those 
which include only the clusters as the 'Single-Plane'. An 
example of the critical curves of a massive cluster with 
and wit hout the line of sight lens p l anes i s shown in Fig- 
ure 3 of IDalal. Hennawi. fc Bodel {2005). There it can 
be seen that the critical curves are only slightly affected 
by the line of sight density fluctuations. Note that in 
order to accurately compute the angular separations of 
the arcs from the cluster, we must take into account the 
deflection of the cluster center, which can be as large as 
~ 1', by lensing from the large scale structure between 
the observer and lens redshift. 

We ray traced through a total of 13,594 unique projec- 
tions of the 878 clusters in our simulation volume for the 
snapshot at Zd = 0.41. The source plane redshift was set 
to z s = 2.0. BCGs were not added to the dark matter 
mass distribution for this analysis. In Figure El we com- 
pare the cumulative distribution of giant arc separations 
for the 'Multi-Plane' and 'Single-Plane' ray tracing sim- 
ulations. It is clear that large scale structure projections 
do not have a significant effect on the abundance of gi- 
ant arcs, as the l argest difference between th e two abun- 
dances is < 7%. IWambsganss et ahl (2004b) reported a 
larger effect, ~ 25%, for a similar source redshift plane, 
based on qualitative arguments regarding the fraction of 
surface mass density pixels which became super-critical 
because of contributions from multiple lens planes. Their 
larger frequency of multi-plane lensing are most likely 
smaller separation, 8 < 10", galaxy scale lensing events 
which would not be considered giant arcs. Our full ray- 
tracing simulations show quantitatively that line of sight 
projections have a modest effect on the statistics of giant 
arcs with separations 9 > 10". 

5. ANALOG HALOS 

To elucidate which properties of CDM halos are most 
important for strong lensing, we have measured struc- 
tural properties for our simulated clusters and generated 
synthetic analog halos which preserve one or more of 
these properties. By comparing the lensing strengths of 
the synthetic analog halos to the N-body halos, we are 
able to quantify the significance of various aspects of halo 
structure for strong lensing. 

The key structural characteristics we measure for our 
halos are the radial density profile p(r) measured with 
respect to the most bound particle, the inertia tensor of 
the particle distribution, and the degree of substructure. 
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We have generated four types of analog halos: 

1. Spherical analogs preserve the radial density pro- 
file p(r) but are spherically symmetric. 

2. Triaxial analogs are similar to spherical analogs, 
however they preserve the inertia tensor of the clus- 
ter core. 

3. Clumpy core analogs are similar to triaxial 
analogs, however the particle distribution within 
the cluster core is left unchanged while the mass 
exterior to the core is smoothed triaxially. 

4. No substructure analogs preserve the particle 
distribution, except for particles within massive 
subhalos, which are triaxially redistributed. 

For the simplest cases of spherical and triaxial 
analogs which we describe below, these halos mir- 
ror the types of simplified analytical models which 
have been used to describe strong lensin g by galaxy 
cluste r s previously in the literature iBartelmannl 
19961 IMeneghetti. Bartelmann. fc Moscardinil 12003a; 
Qguri et all 12003 IQguri fc Keetonl 12004^ It is well 
known that thes e simple models underpredict the number 
of gia nt arcs l)Meneghetti. Bartelmann. fc Moscardinil 
I2003a|) . Our analog halos have the added advantage that 
they exactly reproduce the joint distribution of halo 
shapes and concentrations in the simulation volume. 
Thus, they represent the case of the more sophisticated 
analytical models which attempt to account for the joint 
distributions of halo shapes and concentration mea- 
sured from nume r ical simulations l|Jing fc Sutd 120021 
IQguri et all 120031 IQguri fc Keetonl I2004J) . Although 
analytical models for the dist ribution of substructure in 
dark matter halos exist (e .g. iShethl 120031: iSheth fc Jainl 
120031 IQguri fc Led l2004|) . computing strong lensing 
cross sections for such models is sufficiently cumbersome 
that it has not yet been undertaken. Our ray tracing 
computation for the no-substructure halos allow us 
to quantify the effect of dense substructure on the 
efficiency for producing giant arcs. Because the clumpy 
core analog halos are roughly identical to their parents 
inside of an average overdensity of ~ 10, 000, but are 
smooth ellipsoids outside of this contour, a comparison 
with their parent halos can quantify the degree to 
which the structure of the large scale mass distribution 
changes strong lensing cross sections. We note that 
some of these issues have already been exp lored in 
IMeneghetti. Bartelmann. fc Moscardinil l)2003a[l . There 
the differences between elliptical (triaxial) and real dark 
matter halos was quantified via a multipole analysis in 
Fourier space. Our method of Analog Halos is com- 
plementary, and because we modify these halos in real 
space, we can pinpoint the particular halo properties 
which are important to strong lensing. 

5.1. Generating Analog Halos 

Spherical analog halos are generated by randomly re- 
distributing the angular coordinates of particles while 
holding fixed their radial coordinates. The radial profile 
p(r) of each spherical analog is thus identical to its par- 
ent halo. Note that all radii are measured with respect 
to the most bound particle in the cluster. 



To generate triaxial analog halos we follow a similar 
procedure as for the spherical case. For each particle 
we co mpute the ellipsoidal radius (see e.g. IJing fc Sutd 
2002) 

i?e = V^ + ^ + ^' 

as well as A = (g2<Z3)~ 1//3 , B = qiA, and C = q%A, where 
the coordinates X, Y, Z are measured along the princi- 
pal axes of the core (i.e. the eigenvectors of the inertia 
tensor Q, see eqn. |HJ|. Each particle is then randomly 
mapped onto the surface of the ellipsoid defined by con- 
stant R e - This procedure prese rves the density p rofile 
in ellipsoidal coordinates p(R e ) ijJing fc Sutdl2002j) and 
also maintains the axis ratios and principal axes of the 
core of the original N-body halo. 

The no-substructure analog halos are generated by re- 
distributing all of the particles linked into FOF subhalos 
{b = 0.05) with masses M FO f > 10 115 h" 1 M Q , except 
for those in the largest subhalo, which are left untouched. 
The subhalo particles are redistributed along the ellip- 
soid defined by R e in eqn. (|ll|l using the principal axis 
directions and axis ratios of the core — i.e. using exactly 
the same procedure as for the triaxial analog halos — thus 
also preserving p(R e )- As discussed in H2.3I the FOF al- 
gorithm with linking length b = 0.05 links particles into 
groups enclosed by an average overdensity of ~ 10, 000. 
These analog halos are thus identical to their correspond- 
ing N-body halos except that all substructures contain- 
ing M > 10 115 M Q within a mean overdensity of 
10, 000 have been smoothed out. 

We generate the clumpy core analog halos in exactly 
the same way as the triaxial case except that the particles 
in the most massive subhalo identified by the FOF with 
b = 0.05, i.e. the cluster core, are left untouched. Thus 
the Clumpy Core analog halos are roughly identical to 
their parents inside of a mean overdensity of ~ 10, 000, 
but smooth and nearly identical to the 'Triaxial' analog 
halos outside of this overdensity. 

Figure compares the cluster lens in our simulation 
with the largest lensing cross section to all of its differ- 
ent analog halos. This cluster, taken from a snapshot 
at z d = 0.41, has mass M vir = 8.1 x 10 14 h^ 1 M Q and 
concentration cnfw = 6.9. The axis ratios of the core 
are q^ = 0.46 and q% = 0.40. In Figure 01 we compare 
the radial density profile of this cluster to the profiles 
of the different analog halos. The spherical case is not 
shown because its radial density profile is identical to 
the parent by construction. Note that although we keep 
the ellipsoidal radius R e fixed when redistributing parti- 
cles to create the triaxial, no substructure, and clumpy 
core analog halos, in practice the radial profiles p(r) are 
nevertheless preserved to a high degree of accuracy. 

5.2. Lensing Comparison 

We generated each type of analog halo for all 878 clus- 
ters in the the snapshot at Zd — 0.41. We ray traced 
through every cluster for a source plane at z s = 2.0, av- 
eraging over the same orientations for each halo (see §0J, 
and computed the cumulative distribution of image split- 
tings following eqn. |7|. Since our aim is to understand 
which characteristics of the dark matter distribution are 
most important for strong lensing, BCGs were not added 
to the mass distribution for this analysis. In Figure [S]wc 
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Fig. 3. — Density distributions of four different types of analog halos compared to the original parent. The surface density for three 
projections are shown on the coordinate axis planes to aid visualization. The spatial axes are in units of comoving h Mpc and both 
density (2-d and 3-d) colormaps are logarithmically stretched. Upper left is the the spherical analog, the upper right is the triaxial analog, 
middle left is the clumpy core analog, middle right is the no substructure analog, and the parent halo is depicted at the bottom. The radial 
mass profiles of the analog halos are compared to that of the parent halo in Figure HI 
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Fig. 4. — Comparison of radial density profile p(r) of analog ha- 
los shown in Figure [3] to original parent halo. The y-axis has units 
of h^ 1 Mq/pc 3 . The black points and error bars are the density 
profile and bootstrapped variances of the original parent halo. The 
solid green line is the best fit NFW profile. The inner and outer 
dashed vertical lines indicate the locations of the scale radius and 
virial radius for the best fit NFW profile, respectively. Red squares, 
magenta triangles, and blue upside down triangles, show the den- 
sity profile of triaxial, and clumpy core analog halos, respectively. 
The spherical case is not shown because its radial density profile 
is identical to the parent halo by construction. Although we kept 
the ellipsoidal radius R e fixed when redistributing particles for the 
analog halos shown, the radial profiles are nevertheless preserved 
to a high degree of accuracy. 
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Fig. 5. — Comparison of the number of giant arcs on the entire 
sky for original and analog halos. The black curve is the cumulative 
distributions of giant arcs lensed by the original 'parent' in the 
snapshot at redshift = 0.41. The source plane was at z s = 
2.0, and the density of background galaxies in the reds hift r ange 
[1.75,2.5] was computed according to the discussion in S 17.11 The 
green (long dashed) curve is the distribution of arc separations if 
all clusters in the simulation volume are replaced by their spherical 
analog halos. The red (dot-dashed), blue (dotted) and magenta 
(short dashed) curves show the same quantity but for triaxial, no 
substructure, and clumpy core analog halos, respectively. The data 
plotted above are also shown in Table HI 



compare the number of giant arcs produced by the real 
and analog halos, assuming an area equal to the entire 



sky. The data used to construct this figure is given in 
Table HI 

From the data in the table, we see that spherical ha- 
los underpredict the abundance of giant arcs by a factor 
as large as 50 and triaxial models fall short by a fac- 
tor as much as 60%. These results are consistent with 
the fin dings of iMeneghetti. Bartelmann. &: Moscardinil 
l)2003a|) . Because the distribution of halo concentrations 
and shapes of the parent halos is exactly reproduced 
in our analog halos, we conclude that even the more 
complicated analytical models which convolve analyti- 
cal cross sections wit h the distributions of halo concen- 
tratio ns and shapes l|Qguri et alJl200.lt iQguri fc Keetonl 
120041) , will underpredict the number of cluster lenses by a 
large amount. It is noteworthy that triaxiality increases 
the number of giant arcs by a factor of 4-25 compared 
to the spherically symmetric halos. This is because the 
shallow density cusps p oc r" 1 of CDM halos result in 
an extreme sensitivity to triaxiali ty, as has been empha- 
sized by several authors (DHH ; iDalal fc Keetonl feOOfl 
iBartehnann fc Meneghettill2004|) . 

From a comparison of the real halos to the no sub- 
structure analogs, it is apparent that projections of 
dense substructure onto the small radii probed by strong 
lensing has a relatively modest effect on the lensing 
cross secti ons, increasing the tota l number of arcs by 
~ 5 — 10%. Meneghetti et all l)2000fl came to similar con- 
clusions about changes in lensing cross sections caused by 
baryonic substructure (cluster galaxies) which have small 
masses < 10 12 hr 1 M©. Here we find a similar result for 
dark matter substructure, which can be much more mas- 
sive for the extreme case of a binary or merging cluster. 
We revisit the effect of substructure on strong lensing 
in § El when we statistically quantify the properties of 
cluster lenses. 

Recall that the clumpy core analog halos are identi- 
cal to the original parent halos inside a mean overden- 
sity of 10, 000, but are smooth ellipsoids outside of this 
contour. The close agreement between the clumpy core 
analogs and the parent halos, indicates that halo sub- 
structure and the large scale mass distribution do not 
significantly influence strong lensing cross sections. The 
number of giant arcs produced by CDM halos is 'primar- 
ily determined by the mass distribution within a mean 
overdensity of ~ 10,000. 

The difference between the cross sections for triaxial 
halos and the full N-body halos indicates that some sys- 
tematic departure from ellipsoidal symmetry, on aver- 
age, enhances lensing strength. Some possibilities for 
this include substructure on very small scales, or per- 
haps lopsidedness in the form of boxiness or diskiness 
in the isodensity contours. Naively, we would not expect 
massive substructures to persist in the high-density cores 
of clusters since the dynamical friction timescale is much 
shorter than the Hubble time. We have attempted to 
detect whether larger-scale departures from ellipsoidal 
symmetry enhance lensing cross section by running a 
principal component analysis on the halo density pro- 
files. The results of this analysis were unclear: the prin- 
cipal components we recovered with high signal to noise 
appeared to correspond to varying concentration and el- 
lipticity, with no other significant components apparent. 
This point warrants future study, since it suggests that 
the abundance of giant arcs is sensitive to the detailed 
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TABLE 1 
Analog Halo Comparison 



Type 


N(6 > 10") 


N(6 > 15") 


N(6 > 20") 


N(9 > 25") 


N(6 > 30") 


Real 


274.6 


172.3 


111.2 


74.4 


48.9 


Clumpy Core 


275.1 


172.6 


111.3 


73.8 


48.5 


No Substructure 


264.7 


163.9 


104.7 


69.1 


44.8 


Triaxial 


219.7 


132.7 


81.3 


50.3 


30.4 


Spherical 


60.3 


36.4 


19.4 


8.0 


1.2 



NOTES — Comparison of the number of giant arcs produced by real halos to the number produced by analog halos. The quantity 
N(0 > 10") is the cumulative number of of giant arcs on the entire sky for a source plane at z s = 2.0. The data in this table are plotted in 
Figure 151 
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Fig. 6. — Lensing (dashed) and total (solid) distributions of the 
reduced \ 2 of the NFW fit to the number profile of each halo. The 
median of both the lensing and total distributions, indicated by the 
arrows, are labeled in the upper right . The fact that the median 
value of the total distribution is near unity gives us confidence that 
the NFW profile is providing a good fit to the simulated clusters 
a nd t he bootstrapped ai are accurately describing the errors (see 
§ 12.21 . The similarity between the lensing and total distribution 
indicates that NFW profiles provide just as good a fit to lensing 
clusters as they do to the total population. 



morphology of dark matter halos on small scales. 

6. CHARACTERIZING THE CLUSTER LENS POPULATION 

In this section we address the following question: How 
different is the population of lensing selected clusters 
from the general cluster population? Because cluster 
strong lensing is sensitive to the mass distribution inte- 
rior to the NFW scale radius r < 300 h^ 1 kpc, it provides 
a means to measure the concentration of clusters, espe- 
cially when combined with larger scale mass measure- 
ments like weak lensing llKneib et al.l2003llGavazzi et al.l 
120031 IBroadhurst et alJl2005b|) . It is thus interesting to 
quantify the degree to which these concentrations will be 
systematically higher than the typical cluster in the Uni- 
verse. Besides concentration, other biases might exist. 
Are lensing clusters significantly more or less triaxial? 
Do they have more/less substructure? How likely is the 
major axis of a lensing cluster to be aligned with the line 
of sight? 

We are in a position to answer these statistical ques- 
tions because of the large number of clusters and ori- 



entations that we have ray traced. Specifically, for the 
snapshot at Zd = 0.41 we ray traced through a total 
of 13,594 unique projections of the 878 clusters in our 
V = (320 hr x Mpc) 3 simulation volume. This consti- 
tutes an increase of two orders of magnitude in both 
the number of clusters and number of orientations of 
any previous study. For the remainder of this section, 
we focus our statistical analysis on the clusters in this 
snapshot at Zd = 0.41. Note that BCGs were added to 
the mass distributions of these clusters as described in 
§ 13.11 and thus arc included in the strong lensing cross 
sections. However, we measured the cluster parameters 
(mass, concentration, triaxiality) from the dark matter 
particle distributions alone. 

We define the probability distribution of the cluster pa- 
rameter u (such as mass, concentration, axis ratio, etc.) 



-Plcns(w) 



1 da c s 
J2 a cS du 



(12) 



where Yl a cff is the total effective cross section of all 
the clusters (see eqn. in our simulation volume and 
^f 5 * is the differential distribution of the statistic u, 
which we compute by summing the total cross section 
with u between u and u + Au. In what follows, we use 
Vcs(0 > 15") , the effective cross section for forming giant 
arcs with separations > 15". A comparison of the lens- 
ing distribution Pi ens (u), to the probability distribution 
of the total cluster population P t ot (u) will indicate the 
degree to which the lenses constitute a biased sample. 

In what follows we will compare the NFW parameters 
of cluster lenses to the total cluster population. This 
comparison is clearly only sensible provided that NFW 
profiles provide a good fit to the cluster lenses. In Fig- 
ure we compare the distribution of the reduced \ 2 01 
the NFW profile fits (see eqn. |3J) for the cluster lenses to 
the those for the total population. The similarity of the 
distributions indicates that NFW profiles provide just 
as good a fit to lensing clusters as they do to the total 
population. Furthermore, the fact the reduced x 2 ~ 1 in- 
dicates that NFW profiles provide good fits to the mass 
distributions of our clusters. 

6.1. Mass and Concentration 

In the left panel of Figure [7| we show a scatter plot of 
the orientation-averaged effective cross sections a c s{9 > 
15") against virial mass, for all 878 of the clusters in the 
snapshot at Zd = 0.41. The colors and sizes of the points 
indicate the value of the concentration of each cluster. 
The histogram shows the strong lensing probability av- 
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Fig. 7. — Le/4: Scatter plot of effective strong lensing cross section for forming giant arcs with 8 > 15" against cluster mass, for clusters 
at 2d = 0.41. Each point corresponds to a cluster in our simulation volume and the cross sections have been averaged over all orientations 
(see §|3J. The y-axis has units of arcsec , and the colors and sizes of the points indicate the value of the lens strength parameters K s (see 
§ 12.21 . For M > 10 15 h^ 1 Mq, all clusters are effective strong lenses; whereas at lower masses M < 3 X 10 14 h,- 1 Mq only the most 
concentrated members of the population can multiply image quasars. The histogram shows the mean strong lensing probability averaged 
over mass bins of width 0.1 in log 10 M. Clusters which have zero strong lensing cross section for all orientations are plotted in the lower 
panel, where we have stretched the y-axis so that the points do not all overlap. Right: Comparison of the probability distribution of the 
virial mass, -Picns(l°gio ^vir) °f cluster lenses (dashed red curve) to the mass function, normalized to unity over the mass range considered, 
(solid black curve). These are the same clusters at = 0.41 shown in the left panel. The bin spacing is the same as in the left panel. The 
red histogram fi e ns(l°gio ^0 can be obtained from the scatter plot in the left pan el by summing the values of the cross sections in each 
mass bin and dividing by the total cross section integrated over all bins (see eqn. 1121 . The mean value of log 10 M of the cluster lens sample 
is M vir = 10< lo sio «)ion B = 4.5 x 10 14 h' 1 M Q , indicated by the arrow. 



eraged over mass bins of width dlog 10 M = 0.1. Sub- 
critical clusters which had zero effective cross sections 
are plotted in the lower panel, where we have arbitrarily 
stretched the vertical axes so the points do not overlap. 
The scatter plot illustrates the range of cross sections 
which exist in a given mass bin. In addition, it shows 
that for the largest masses characteristic of superclusters, 
Al ~ 10 15 h^ 1 M Q , all clusters are effective strong lenses, 
whereas at lower masses M < 3 x 10 14 h^ 1 M Q , only 
the more concentrated members of the population are 
strong lenses. Notice the extremely steep dependence of 
the lensing cross section on mass: the mean cross section 
changes by nearly four orders of magnitude over a single 
decade in mass. In the right panel, we compare the prob- 
ability distribution Pi ens (log 10 M) (dashed red curve) to 
the cluster mass function normalized to unity over the 
mass range M > 10 14 h^ 1 M Q considered (solid black 
curve). The bin spacing in the right panel is the same 
as in the left panel. The red histogram Pi cns (log 10 M) 
can be obtained from the scatter plot in the left panel by 
summing the values of the points in each mass bin and 
normalizing the histogram by the total number of lenses 
(see ean. 112(1 . The median value of log 10 M for the cluster 
lens sample is M = 10 lo Sio^„ 5 = 4 5 x 10 14 ft-i m q . 

The left panel of Figure0suggests that lensing selected 
cluster samples are likely to show significant concentra- 
tion bias. Since we have seen that cross section is such 
a steep function of cluster mass, it is desirable to first 
normalize out the mass dependence of cluster concen- 



tration. Thus we consider the distribution of the ratio 
c v i r /c v ir(Af), where c v i r (M) is some kind of average con- 
centration. Because simply computing the mean concen- 
tration in mass bins will be sensitive to outliers and very 
noisy at the high mass end where we have few clusters, we 
must be careful about the fitting procedure. In Figure 03 
we show a scatter plot of concentration versus virial mass 
measured from the clusters in our simulation volume at 
Zd = 0.41. The blue curve is a second order B-splinc 
fit to points with breakpoints set every 175 clusters or- 
dered by mass. This fit was computed iteratively with 
3cr outliers rej ected at each i terati on until convergence 
was achieved. Bullock et al.l l|2001|) measured the trend 
c = (!^ z ) (M/M^)~°- 13 from a large ensemble of numeri- 
cally simulated clusters, where M* is the nonlinear mass 
at z = which is M* = 1.3 x 10 13 h" 1 M Q for our cos- 
mology. The red dashed curve in Figure El c urve shows 
the result of adopting the lBullock et alJ l(2001f ) power law 
scaling and fitting the median concentration in five log- 
arithmically spaced mass bins for the linear amplitude. 
We get c = (1^) (M/M»)~°- 13 , i.e. our average concen- 
trations are 37% higher than lBullock et all l)2001|) 7 . Our 

7 This disagreement is likely due to differences in the halo fitting 
procedure and also possibl y because of the different cosmological 
model simulated. IBullock et al.l 120011) truncate halo radii when 
halos overl ap a nd they only fit particles bound to the halo. We 
argued in S I2.2l that neither practice is appropriate to the context 
of gravitational lensing. 
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B- spline fit is in near a greement with the relation found 
bv lBullock et al.1 (|2001f ) - although the mass scaling is not 
quite as steep. 

The lensing and total distributions of c/c(M) are 
shown in the right panel of Figure |H1 The median of the 
total distribution is consistent with unity to better than 
1%, indicating that the average concentration shown in 
the left panel of Figure [HI is a good fit to c(M) . We see 
that the cluster lens population has three dimensional 
concentrations ~ 18% higher than the average cluster at 
a similar mass. 

6.2. Triaxility and Orientation Bias 

In this section we investigate alignment bias and study 
the effect of triaxiality on strong lensing probability. The 
lensing and total distributions of the absolute value of the 
cosine of the angle between the line of sight and the major 
axis direction of the cluster, |cos#| = z • V a , are shown 
in the left panel of Figure The thin (black) line is 
the distribution for all clusters, which is flat as expected. 
The axis ratio qi is divided into three groups with equal 
total lensing probability. The magenta (dotted), blue 
(dashed), and green (dot-dashed) curves are individual 
lensing contributions from clusters with axis ratio q2 in 
the lower third q2 < 0.50 (most triaxial), middle third 
0.50 < (72 < 0.66, and upper third qi > 0.66 (least tri- 
axial). These contributions sum to give the total lensing 
distribution, which is the thick red (solid) curve. The 
median value for the lensing clusters is | cos#|i on s = 0.67, 
indicating significant alignment bias. Furthermore, the 
alignment bias is larger for the more triaxial clusters, as 
expected. Our findings disagree with t he discussion in 
iBartelmann. Steinmetz. fc Weiss! l|1995j) who also com- 
pared the orientation of the cluster principal axes to the 
line of sight and found no significant correlation between 
the two. Because this study only ray traced 3 different 
projections through 13 clusters, it is likely that they did 
not have the statistics to measure the correlation which 
we detect with high significance. 

Despite this large tendency for lensing clusters to be 
aligned with the line of sight, the lensing population has 
nearly the same distribution of axis ratios as the total 
cluster population. This is the clear indication of Fig- 
ure which shows the distributions of the axis ratios 
q2 and ^3 (shown in the left and right panels, respec- 
tively). The lensing distributions are nearly indistin- 
guishable from the total cluster population. 

We have detected significant alignment bias between 
the line of sight and the principal axis direction of the 
clusters, and furthermore we argued in [|S]that triaxial- 
ity significantly increased the total lensing cross section 
for CDM clusters. On the other hand, Figure indi- 
cates that lensing clusters have the same distribution of 
axis ratios as the total cluster population. These two 
statements seem to contradict each other. If triaxiality 
enhances lensing cross sections and if their is a tendency 
for strong lenses to be aligned with the line of sight, why 
are lensing clusters not systematically more triaxial than 
the total population? 

One possible explanation for this apparent contradic- 
tion is that the shape of a dark halo is correlated with 
its concen tration, as mi g ht be expected on dynamical 
grounds. lAllgood et alJ l)2005|) have found that simu- 
lated dark halos which have recently collapsed are more 



triaxial ((72 smaller) and less concentrated, whereas ha- 
los which have collapsed in the more distant past are on 
average less triaxial (92 larger) and more concentrated. 
Dynamically old clusters formed when the Universe was 
more dense and hence are more concentrated, and sim- 
ilarly have had more time to relax and hence should be 
less triaxial. Indeed, a similar trend exists in our sample 
of clusters as is shown in the scatter plot of normalized 
concentration versus axis ratio in Figure ITD 

6.3. 2-D vs. 3-D concentration 

Because strong lensing probes the mass in projection, 
the three dimensional concentration considered above is 
not an observable. Rather, the two-dimensional concen- 
tration is measured from detailed modeling of the arcs 
and image positions. For a given 3-D concentration, we 
expect an additional bias to exist in the distribution of 
2-D concentrations since lensing clusters are triaxial and, 
as discussed above, tend to be viewed down their major 
axes. The left panel of Figure ^| indicates that this is 
indeed the case. There we show the lensing and total 
distributions of C2d/czd, which is the two-dimensional 
concentration (for a given projection) normalized by the 
three-dimensional concentration of each cluster. The his- 
togram is computed using the cross section for each of 
the 13,594 unique projections of the 878 clusters in our 
simulation volume. Note that median value of the total 
distribution deviates from unity by 6%. This bias rel- 
ative to the three dimensional concentrations is caused 
by the fact that our two dimensional NFW profile fits 
underestimate the concentrations because we fit to the 
projection of a profile which extends to infinity (see dis- 
cussion in § 12. 2J) . Nevertheless, we see that the lensing 
distribution of C2d/c3D is biased by ~ 19% relative to 
the total distribution. 

We compare the lensing and total distributions of the 
two dimensional concentrations in the right panel of Fig- 
ure El The concentrations of lensing clusters are 34% 
larger than than those of the total cluster population. 
This bias is caused by the combination of two effects. 
First, strong lensing clusters have three dimensional con- 
centrations which are on average 18% larger than the 
typical cluster at the same mass, as is indicated by Fig- 
ure [S] Second, given a three dimensional concentration, 
Figure IT21 indicates that strong lensing prefers the pro- 
jections through clusters which give ~ 19% higher con- 
centrations. This 34% concentration bias should be kept 
in mind when comparisons are made between concentra- 
tions measured from modeling observed strong lenses to 
the mean concentrations measured from N-body simula- 
tions. 

The foregoing discussion has bearing on the recent 
measurements of anomalously high concentrations from 
detailed modeling of individual lensing clu sters. From 
combi ned weak and strong lensing analyses. iKneib et alJ 
(2003) measured a best fit concentrati o n of c nfw = 22 
for CL0024+1654, while lOavazzi et all l|2003f) measured 
cnfw = 12 for MS 2137.3-2353 (but s ee lDalal fc Keetonl 
2003), and iBroadhurst et al.l (|2005bf) measured c v - 1T = 
13.7 for Abell 1689. It is unlikely that the brightest clus- 
ter galaxies in these clusters increase the concentrations 
over our expectation for dark matter alone. 8 Note that 

8 For CL0024+1654 and Abell 1689 the arcs occur at large radii 
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Fig. 8. — Left: Scatter plot of cluster concentration versus virial mass for all the clusters in our simulation volume at = 0.41. The 
solid blue curve is a second order B-spline fit to the points with breakpoints set every 175 clusters ordered by mass. To reduce sensitivity 
to outliers, this fit was computed iteratively with 3<r outliers rejected at each iteration until convergence was achieved. The dashed red 
curve shows the fit, c v i r = T^p-fM/M,) -0 ' 13 , which is the result of adopting the Bullock et al. (2001) power law scaling and fitting the 

median concentration in five log arithmicall y sp aced mass bins for the amplitude. The B-spline fit (solid blue line) is in near agreement 
with the relation found by Bull ock et alj 12001). although the mass scaling is not quite as steep. Right: Lensing (dashed) and total (solid) 
distributions of the ratio c/e(M). Arrows indicate the median of the distributions. The fact that the median of the total population is 
consistent with unity indicates that our fit for c(M) has effectively removed the average mass scaling from the concentration. Cluster lens 
have three dimensional concentrations 18% higher than the typical cluster with similar mass. 




0.0 0.2 0.4 0.6 0.8 1.0 

|COS0| 

Fig. 9. — The lensing and total distributions of the absolute value 
of the cosine of the angle between the line of sight and the major 
axis direction of the cluster. The thin (black) line is the intrinsic 
distribution, which is flat as expected. The axis ratio <j2 is divided 
into three groups with equal lensing probabilities. The magenta 
(dotted), blue (dashed), and green (dot-dashed) curves are indi- 
vidual contributions from clusters with axis ratio <j2 in the lower 
third Q2 < 0.50 (most triaxial), middle third 0.50 < f?2 < 0.66, and 
upper third 52 > 0.66 (least triaxial). These contributions sum 
to give the total distribution which is the thick solid (red) curve. 
The median of both the lensing and total distributions, indicated 
by the arrows, are labeled on the plot. The median for the lensing 
clusters is |cos0| = 0.67, indicating significant alignment bias. 



> 30", where the baryonic component contributes only a small 
fraction < 10% of the total mas s enclosed by the critical curves 
(see e.e. fBroadhurst et alj|2005al) . Baryons are a larger concern 
for MS 2137.3—2353 because its arcs are at smaller radii < 15", 
however Gavazzi et al. (2003) found that including a mass compo- 



iKneib et al] (|2003j) and lGavazzi et alJ l|2003j) use a con- 
vention for the NFW concentration, cnfw = r 20o/ r s, 
where r2oo is the radius at which the average density 
is 200 times the critical density, which differs from our 
definition of c v ir = r V i r /r s , where r V i r is the radius 
where the average density is A v i r (z) times the mean den- 
sity. Converting these concentrat ions to our convention 
<|Whitei200HlHu fc Kravtsovl2003j) . gives c vir = 26.2 and 
Cvir = 14.6, for CL0024+1654 and MS 2137.3-2353, re- 
spectively. 

These high concentrations are very puzzling, consider- 
ing that the distribution in the right panel of Figure 1121 
predicts that the probability for C2D > 14 is less than 
2%. Why should the three best studied lensing clus- 
ters in the Universe all have concentrations on the tail 
of the concentration distribution? One possible expla- 
nation is that these clusters have been studied inten- 
sively because they are known to contain multiple gi- 
ant arcs, and one might expect clusters with multiple 
high-surface brightness arcs to be even more biased than 
the gener al population of clus ter lenses. While not im- 
plausible, IHo fc White] <|2004|) found that ~ 40 - 50% 
of massive clusters that produce giant arcs show multi- 
ple arcs, which suggests that multi-arc clusters are rel- 
atively common. A definitive solution to this puzzle of 
high observed concentrations will clearly require a large, 
homogeneously selected statistical sample of strong lens- 
ing clusters. 

6.4. Substructure 

In § 12.31 we defined three simple statistics which quan- 
tify the amount of substructure in a cluster. The first, 

nent associated with the central galaxy did not significantly change 
the resulting concentration. 
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Fig. 10. — Lensing (dashed) and total (solid) distributions of the axis ratios. The left panel shows the distributions of the axis ratio 
52 and the right panel shows the distributions of 53. The mean of both the lensing and total distributions, indicated by the arrows, are 
labeled in the upper right corner of each plot. The lensing population has nearly the same distribution of axis ratios as the total cluster 
population. 
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Fig. 11. — Scatter plot of normalized concentration, c/c(M) 
versus axis ratio 52 for the clusters in our simulation at = 0.41. 
The black triangles show clusters with zero strong lensing cross 
section <T cf f(> 15") and circles show clusters with nonzero cross 
sections, where the colors and sizes of the points indicate the size 
of the cross section in units of arcsec 2 . Cluster with a large degree 
of triaxiality (small 92) tend to be under concentrated; whereas, 
clusters which are less triaxial (large 52) tend to be over concen- 
trated. This correlation between 92 and normalized concentration 
is expected on dynamical grounds (see discussion in Sl6l. 



-^sub-core, is the ratio of the mass in substructures to 
the core mass of the cluster; the second, iVf B ub— vin is the 
ratio of the mass in substructure to the virial mass of the 
cluster; and the third, M2- V i r , is the ratio of the most 



massive substructure (i.e. second to the the cluster core) 
to the virial mass. The lensing and total distributions of 
these three statistics are shown in Figure IT31 

A comparison of the lensing and total distributions of 
these statistics indicates that any correlation between 
strong lensing cross section and the presence of substruc- 
ture is marginal at best. First, focusing attention on 
Figure^D which is for substructures within the virial ra- 
dius, the median of the statistic Af su b_ C orc for the lensing 
and total distributions suggests that lensing clusters have 
~ 10% less substructure than the total cluster popula- 
tion. The mass in substructure normalized by the virial 
mass, M su b-vir, is biased in the opposite direction by 
~ 10% compared to the total cluster population. Both of 
these statistics consider all substructures > 10 12 h -1 M Q 
(see eqn. |SJ) and effectively integrate over the subhalo 
mass function. In contrast, the statistic M su b-2 quan- 
tifies the binarity of the cluster and is expected to be 
largest for a cluster undergoing a major merger. The 
median value of M su b-2 for lensing selected clusters is 
indistinguishable from the total cluster population. 

The fact that M su b-2 does not correlate with strong 
lensing brings into question the notion that unrelaxed 
clusters undergoing major mergers are more e f fective 
gravitational lenses. In particular, iTorri et al.l l)2004f) 
analyzed two different projections through a single tem- 
porally resolved cluster merger and argued that mergers 
significantly enhance strong lensing cross sections. While 
individual projections can be enhanced by the presence 
of significant substructure and similar examples exist for 
our clusters, these enhancements due to chance align- 
ments of substructures are diluted, in the orientation av- 
erage, by the other projectio ns for which substructures 
are not aligned. In addition, iZentner fc Bullockl l)2003f) 
found that the degree of substructure in dark matter ha- 
los anti-correlates with concentration, and we have seen 
that more centrally concentrated clusters tend to dom- 
inate the lensing cross section. Finally, our conclusion 
here that substructure does not correlate with strong 



CLUSTER LENSES 



15 



( -'2D^3D — 0-94 
( -'2D^3D = 1 -1 1 



0.5 




2.0 2.5 3.0 



^2D^3D 




Fig. 12. — Left: Lensing (dashed) and total (solid) distributions of the ratio C2D/C3D where the two dimensional concentrations have been 
measured for each projection through the clusters. Arrows indicate the median of each distribution. Strong lensing prefers the projections 
through clusters which give ~ 19% higher concentrations. Right: Lensing (dashed) and total (solid) distributions of the two dimensional 
concentrations, C2D> measured from each projection through the clusters. Arrows indicate the median of the distributions. Cluster lenses 
have two dimensional concentrations which are 34% higher than the typical cluster in the Universe. 
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Fig. 13. — Lensing (dashed) and total (solid) distributions of the substructure statistics. The left panel shows the distributions of 
Af su b_ corc , which is the mass in substructure normalized by the dense core mass. The middle panel shows M su t,— v in which is the mass in 
substructure normalized by the virial mass. The right panel shows M su b_2, which is ratio of the the single most massive substructure to 
the virial mass. The median of both the lensing and total distributions, indicated by the arrows, are labeled on each plot. The agreement 
between the lensing and total distributions for these three different substructure statistics clearly indicates that the population of cluster 
lenses are no more relaxed or disturbed than typical clusters in the Universe. 



lensing is reinforced by the results of the no substruc- 
ture analog halos in § There we saw that smoothly re- 
distributing all the mass in substructure had a marginal 
effect, ~ 5 — 10%, on the total number of giant arcs 
produced. Our statistical analysis of 13,594 different ori- 
entations of 878 clusters clearly indicates that the popu- 
lation of cluster lenses are no more relaxed or disturbed 
than typical clusters in the Universe. 

7. DO THE ARCS DISCOVERED IN THE RCS AGREE 
WITH ACDM? 

In this section we revisit the issue of the abundance 
of giant arcs in the Re d Cluster Sequence cluster survey. 
Clad ders et alJ l|2003j) surveyed 90 deg 2 and found sev- 
eral giant arcs, all in clusters at high redshift z > 0.6. 
DHH were unable to determine whether the number of 
giant arcs discovered in the RCS is consistent with the 
ACDM model because the si mulation cube used in tha t 
study was too small (but see iWa mbsgan ss" et al.l l2004a). 
Here we revisit this comparison using a simulation vol- 



ume (320/141. 3) 3 = 11.6 times larger. In Figure ITU we 
plot the ratio of the RCS volume to the volume of our 
320 h" 1 Mpc simulation cube as a function of redshift. 
Our simulation volume is larger than the RCS for z ~ 0.7, 
where three of the five RCS arcs which we consider oc- 
cur (see below), and comparable to the RCS volume for 
higher redshifts 0.7 < z < 1.0 where the other two RCS 
arcs land. Hence, we have a sufficient volume to sample 
the very rare massive cl usters likely to be re sponsible for 
the RCS strong lensing l|Ho & Whiteil2004D . 

7.1. The Number Density of Background Galaxies 

Gravitational lensing conserves surface brightness but 
does not conserve total integrated flux. Because arcs are 
resolved, at least in the tangential direction, arc detec- 
tion should be limited by surface brightness rather than 
integrated flux. In the regime where arcs are resolved 
both tangentially and radially, there is no magnification 
bias. However, high redshift galaxies are known to be 
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Fig. 14.— The ratio of volume of the 90 deg 2 RCS survey 
to our simulation volume as a function of redshift. The quantity 
(Vj/L 3 ) is plotted for each snapshot redshift Zj where Vj is given 
by eqn. an d L = 320 /i -1 Mpc is the size of our simulation cube. 
The snapshot at redshift Zj, was taken to represent the volume of 

the Universe over the redshift range 1 - 2 -^- s ^\, which is 

represented by the bins of the histogram. 
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Density of Background Sources 
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NOTES — Density of background sources used in this paper. 



comp act with half light radii r h ~ 0.3 ({Ferguson et al.1 
|2004[) and thus giant arcs could be radially unresolved. 
This depends on both the radial magnification of the 
lensed image and the res olution or seeing. Ab out half 
of the arcs published by iGladders et alJ 1)2003(1 are re- 
solved radially, while the other half are not. For these 
radially unresolved arcs there is a small amount of mag- 
nification bias because the image is integrated over the 
seeing disk in the radial direction. Because the typical 
radial magnifications will be small (~ seeing/r^ ~ 2), we 
chose to neglect this effect here as it will be negligible 
compared to our dominant uncertainty which turns out 
to be the counts of background galaxies. A more careful 
analysis could include this radial magnification bias by 
assuming a typical seeing and including a distribution of 
source sizes in the Mo nte Carlo part of the ray tracing 
simulation (see e. g.lHo fc Whitell200l . 

Gla dders et al.l l|2003') impose a surface brightness cut 
on their arc sample of \ib_ c < 24. Thus, in order to com- 
pare statistics, we require the surface brightness func- 
tion of high redshift galaxies, rather than the luminos- 
ity function. Note that because the size distribution of 
galaxies is broad ( Fer guson et al.l f2004|) and the lumi- 
nosities and sizes of high redshift galaxies are correlated 
ijBouwens et alJl2004 iTruiillo et alJl2004(l . simply using 
the mean size to convert a limiting surface brightness to 
a limiting magnitude will not give the correct number 



c ounts. 

iFontana et alJ l)2000D publish a photometric redshift 
catalog which inclu des the N e w Te chnology Telescope 
(NTT) Deep Field. IPoli et alJ {l999) measured the half 
light radii of all the galaxies in the NTT deep field using 
image deconvolution techniques. The area of the NTT 
deep field is only 4.8 arcmin 2 . Even for the relatively 
coarse redshift bins we use here Az ~ 1.0, cosmic vari- 
ance for such a small field can be as large as ~ 50% 
ijSomerville et alJ 12004) . Furthermore, half light radii 
measured from space based imaging would be clearly 
preferable to deconvolved ground based data. However, 
no wide field space based catalog which includes redshifts 
and gal axy sizes ex i sts at the time of wri ting. Thus, we 
use th e IPoli et ahl (|1909f) sizes with the IFontana et alJ 
( 2000) photo-z's to determine the number of background 
galaxies in each redshift bin with hr < 24. The density 
of background sources used for each source redshift bin 
are listed in Table [3 We caution that uncertainty in the 
density of sources is likely to be a significant source of 
error in our comparison of the abundance of giant arcs 
to theory. For example, the NTT deep field contains 23 
galaxies in the redshift range [2.5,3.5] with R < 25, or 
4.8 galaxies per arcmin -2 . Yet the incompleteness cor- 
rected luminosity functi on derived from the z ~ 3 Lyman 
break galaxy sample of Stci del et al.l (|1999f) predicts 1.9 
galaxies per arcmin -2 in this range. 

7.2. Comparison with RCS 

The cumulative distribution of image splittings pre- 
dicted for the RCS survey from our ray tracing sim- 
ulations with BCGs added is shown in the left panel 
of Figure 1151 Contributions from the individual source 
planes are also shown, along with the total num- 
ber predicted from dark matter alone. BCGs sig- 
nificantly increase the total number of arcs for sep- 
arations 0__<_^lj)^__consistaxt_w^^ results 
(DHH: IMeneehetti. Bartelmann. fc Moscardml 1200311 
IHo fc Whitell2004D 

We restrict our comparison to giant arcs with separa- 
tions 9 > 10", since for smaller separations, our results 
will be very sensitive to the details of how we paint BCGs 
onto dark halos. We also restrict attention to arcs with 
length to width ratios L/W > 10 . Five of the arcs from 
Table 1 of IGladders et ail l|2003|) satisfy these criteria. 
These are the two arcs in RCS 0224.5-0002 (z = 0.77), 
one arc in RCS 1324.5+28245 (z = 0.85), one arc in 
RCS 1419.2+5326 (z = 0.64), and one of the arcs in 
the secondary high redshift sample RCS 2319.9+0038 
(z = 0.91; Gladders private communication January 
2005). 

The right panel of Figure ^| compares the predicted 
redshift histogram of arcs with 9 > 10" to the redshift 
histogram of the RCS lensing clusters using the same 
binning. It is rather surprising that the RCS did not 
detect any arcs with z < 0.63, as the simulations predict 
~ 5 arcs in this redshift range. For redshifts z > 0.63, the 
simulations predict ~ 4 arcs, whereas all five of the RCS 
arcs are in this range. Although it is odd that the RCS 
arcs all pile up between z = 0.6 — 0.9, this anomaly is 
not statistically significant considering the small numbers 
of objects. We conclude that, at present, there is no 
significant excess of high redshift lensing clusters in the 
RCS survey. 



CLUSTER LENSES 



17 



A 




15 20 25 
[arcsecs] 



3.0 

2.5 

5 — 2.0 
o 

a. 1.5 

1.0 

0.5 
0.0 



| RCS 
| Total 

: z s =1.0 
! z s =1 .5 

! z s =2.0 
I z s =3-0 
I z o =4.0 



0.0 



0.5 



1.0 



1.5 



Fig. 15. — Left: Predicted cumulative distribution of giant arc separations for the RCS cluster survey. The green (dotted), blue (short- 
dashed), orange (dot-dashed), magenta (dot-dot-dashed), and cyan (long dashed) curves are the individual contributions from the source 
planes at z B = 1.0, 1.5, 2.0, 3.0, and 4.0, respectively. The sum of these curves gives the total number of multiply imaged quasars which is 
the thick black (solid) curve. The thin solid (black) line shows the total number of lenses if BCGs are neglected and we ray trace through 
dark matter only. Right: Redshift histogram of cluster lenses with splittings 9 > 10" in the RCS survey. The individual contributions from 
each source plane are also plotted with the same line and color scheme as in the left panel. The thick black (solid) histogram is the total 
num ber of giant arcs predic ted including BCGs. The red (solid) histogram is the observed redshift distribution of cluster lenses in the RCS 
from lGladders et alj J2003). 



Before we conclude this section, we compare the pre- 
dicted number of arcs here to the previous results in 
DHH. They predicted a total of ~ 3 arcs in the RCS 
for z < 1.0, whereas here we predict ~ 8 arcs in this 
redshift range. First, the cosmological parameters for 
the simulations used in that study were slightly differ- 
ent, most importantly the simulation used by DHH had 
er 8 = 0.9, whereas that used in this study used as = 0.95. 
This changes the abundance N(> M) of massive clusters 
by ~ 35 - 40% for M = 3 x 10 14 hr 1 M Q and as much 
as - 170% for M = 3 x 10 15 h^ 1 M Q at z ~ 0.5. Al- 
though significant, the change in cluster abundance does 
not account for the entire disagreement. Another dif- 
ference which could be important for the higher redshift 
(z > 0.6) clusters is the additional high redshift source 
planes (z s = 3.0 and z s = 4.0) used here, whereas the 
highest redshift source plane considered by DHH was 
z s = 2.0. However, it is likely that the the factor of 
11.6 smaller volume used in that study is responsible for 
most of the discrepancy. As indicated by Figure the 
mass function is a very steep function of cluster mass 
and furthermore the scatter about this mean relation is 
large because of the large underlying scatter in cluster 
properties (i.e. Figure |HJ). Large cosmological volumes 
are thus required before convergence to the cosmic mean 
can be achieved. 

Have the results of this study converged to the cor- 
rect cosmic mean lensing cross sections? First note that 
the small volume simulated is at least partly compen- 
sated by averaging over a large number of orientations for 
each cluster. Although the lensing distribution of cluster 
masses in the right panel of Figure [7| indicates that halos 
M ~ 5 x 10 14 h^ 1 Mq are dominating the total lensing 
cross section, the tail of this distribution suggests that 
we may not yet have converged at the high mass end. 
Possible evidence for a lack of convergence is the fact 



that iLi et al.l l)2005|) find much lower optical depths for a 
simulation volume of similar size. Besides convergence to 
the mean, another serious issue is cosmic variance. Since 
the volume we simulated is comparable to the volume 
of the RCS, how likely are these two realizations to give 
the same answer? Answering this question requires one 
to measure the highly non-Gaussian probability distribu- 
tions of very rare events. This is clearly beyond the scope 
of the present study, yet it is a question which must be 
tackled if giant arc statistics are to become a quantitative 
tool for cosmology. 

8. SUMMARY AND CONCLUSIONS 

We attempted to isolate which properties of CDM clus- 
ters make them effective gravitational lenses, by intro- 
ducing several different types of 'Analog Halos,' which 
retain one or more of the properties of the real simulated 
clusters. The results of this analysis are: 

- Spherical halos underpredict the abundance of gi- 
ant arcs by a factor as large as 50 and triaxial mod- 
els fall short by a factor as much as 60% (see Ta- 
ble HJ. 

Triaxiality increases the number of giant arcs by 
a factor of 4-25 compared to the spherically sym- 
metric halos, because the shallow density cusps 
p a r" 1 of CDM halos res ult in an extreme sen- 
sitivity to triaxiality (DHH; iDalal fc KeetorJl2003l 
iBartelmann fc Meneghettill2004D . 

- Projections of halo substructure onto small radii 
and the large scale mass distribution of clusters do 
not significantly influence strong lensing cross sec- 
tions. The number of giant arcs produced by CDM 
halos is primarily determined by the mass distribu- 
tion with a mean overdensity of ~ 10, 000. 
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- The clumpy cores ol dark matter halos result in 
~ 25 — 60% more giant arcs than smooth ellipsoids, 
which suggests that the abundance of giant arcs can 
probe the small scale morphology of dark matter 
halos. 

We measured the properties of a large ensemble of 
clusters and characterized the cluster lens population by 
computing cross sections for each cluster with ray tracing 
simulations. Our statistical comparison of the properties 
of the lensing and total cluster populations yielded the 
following results: 

- NFW profiles provide just as good a fit to lensing 
clusters as they do to the total population of clus- 
ters. 



The typical mass of a lensing cluster is M v - a 
10 14 h~ l Mr. 



4.5x 



- Lensing clusters have 34% higher concentrations 
than the typical cluster at the same redshift with a 
similar mass. This bias is result of a combination of 
two effects. First, the lensing population is biased 
towards clusters with higher three dimensional con- 
centrations. Second, given a three dimensional con- 
centration, orientation bias will favor projections 
along the major axis with higher two dimensional 
concentrations . 

- The anomalously high concentration s c > 14 re- 
cently reported by several groups llKneib et aD 
12003 iGavazzi et al J200l iBroadhurst et alJl2 005b) 
appear inconsistent with the concentration distri- 
bution in our simulations, which predict that < 2% 
of lensing clusters should have concentrations this 
high. 

- The population of cluster lenses are no more re- 
laxed or disturbed than typical clusters in the Uni- 



- Strong lensing clusters tend to have their principal 
axis aligned with the the line of sight. The median 
angle is |cos#| = 0.67. 

- The distribution of axis ratios of strong lensing 
clusters is indistinguishable from the total cluster 
population. 

We revisited the question of whether there is an ex- 
cess of giant arc s detected for high r edshift clusters in 
the RCS survey ijGladders et al.ll2003^) . Our simulations 
predict 9 total arcs in the RCS survey, whereas five were 



discovered. At low redshift (z < 0.6) the RCS found zero 
arcs, which is discrepant with our prediction of ~ 5. At 
high redshift (z > 0.6), our prediction of ~ 4 giant arcs 
is consistent with the five discovered in the RCS. There 
is no significant excess of high redshift lensing clusters in 
the RCS survey over predictions from the ACDM model. 

Finally, we emphasize that the results in this work 
were all based on dissipationless dark matter only N- 
body simulations. Thus, we implicitly assumed that 
dissipative baryonic processes such as heating, cooling, 
turbulence, and star formation have a negligible effect 
on the average surface mass density of a cluster within 
^10,000 ~ 300 hr x kpc (although we di d ac count for the 
effect of brightest cluster galaxies in § 13 - 1 1> . While hy- 
drodynamical simu lations of cluster formation are mak- 
ing rapid progress llGnedin et all2004l:lKazantzidis et alJ 
120041 iKravtsov et alJ l2005|) . thev are not yet mature 
enough to make definitive predictions about the effects 
on the mass distribution on scales relevant for strong 
lensing (Puch wein et al.l|2005t) . This is a fruitful topic 
for future research, considering the potentially significant 
impact these processes could have on our interpretation 
of strong lensing in clusters. 
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